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Abstract 

In some supersymmetric models, the gluino (g) is predicted to be light and 
stable. In that case, it would hadronize to form R-hadrons. In these models, 
the missing energy signature of the lightest supersymmetric particle is no longer 
valid, even if R-parity is conserved. Therefore, such a gluino is not constrained 
by hadron collider results, which looked for the decay g — > qqx?. 
Data collected by the DELPHI detector in 1994 at 91.2 GeV have been analysed 
to search for qqgg events. No deviation from Standard Model predictions is 
observed and a gluino mass between 2 and 18 GeV/c 2 is excluded at the 95% 
confidence level in these models. Then, R-hadrons produced in the squark 
decays were searched for in the data collected by DELPHI at the centre-of-mass 
energies of 189 to 208 GeV, corresponding to an overall integrated luminosity 
of 609 pb _1 . The observed number of events is in agreement with the Standard 
Model predictions. Limits at 95% confidence level are derived on the squark 
masses from the excluded regions in the plane (mq,,mg): 

> 90 GeV/c 2 , and mg i > 96 GeV/c 2 for purely left squarks. 

> 87 GeV/c 2 , and > 82 GeV/c 2 independent of the mixing angle. 



(Eur. Phys. J. C26 (2003) 505) 
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1 Introduction 

In minimal supergravity supersymmetry models (mSUGRA), the gaugino masses (Mj) 
are usually supposed to evolve from a common value at the GUT scale. In such 
models, the Mj are proportional to the corresponding coupling constants (gi) and the 
gluino is naturally heavier than the other gauginos at the electroweak scale. 

Mi M 2 M 3 

-4- = —t = -4- that is, M l : M 2 : M 3 ~ 1 : 2 : 7 

91 9 2 9i 

Nevertheless, models exist where the Mj do not follow this relation. M3 could be lighter 
than the other gaugino masses [1], and in this case, the gluino is the Lightest Supersym- 
metric Particle (LSP). For example there is a particular Gauge Mediated Supersymmetry 
Breaking model (GMSB) [2] where the gluino can either be the LSP or the next to light- 
est supersymmetric particle (NLSP) with a gravitino LSP. In the latter case, the lifetime 
of the gluino would ever be sufficiently large to consider the gluino as a stable particle 
for collider physic. If R-parity is assumed, the gluino is stable in all these models and it 
should hadronize to form R-hadrons because of color confinement. 

The gluino has been intensively searched for in hadron collisions in various decay chan- 
nels [3]. However, the limit obtained [m- g > 173 GeV/c 2 for rrig = m 9 -) does not apply 
to a stable gluino. For this model, it has been shown that CDF run I data could not 
constrain a stable gluino with mass lower than 35 GeV/c 2 [1,2]. On the other hand, the 
gluino mass could be much larger than it is in the so-called light gluino scenario [4] , which 
seems to be excluded by the measurement of the triple gluon coupling and of the four-jet 
rates at LEP [5]. 

A pair of gluinos can be produced in the splitting of a gluon. Figure 1 shows the 
Feynman diagram of this process and the corresponding cross-sections at centre-of-mass 
energies of 91.2 GeV (LEP1) and of 200 GeV (LEP2). The production rate is too low at 
LEP2, and LEP1 data must be analysed to be sensitive to this process. In this channel, 
the gluino does not originate from the decay of another sparticle, so it can be produced 
even if the other supersymmetric particles are not accessible. 

Gluinos can also be produced from other sparticle decays. cjx and are the su- 
persymmetric partners of the left-handed and right handed quarks. With large Yukawa 
coupling running and important off-diagonal terms the supersymmetric partners of top 
and bottom quarks are the most probable candidates for the charged lightest supersym- 
metric particle. The squark mass eigenstates are parametrised by a mixing angle 0q. The 
lighter squark is given by qi = q^sin^ + q^cos^. The stop (ti) and the sbottom (bi) 
could be pair-produced at LEP via e + e~ annihilation into Z°/7. A squark mixing angle 
equal to zero leads to the maximal squark pair production cross-section, while the mini- 
mal cross-sections are obtained for a mixing angle of 56° for the stop and of 68° for the 
sbottom. For these particular angles, the Z° — qi — q\ coupling is suppressed. This paper 
also describes the search for R-hadrons from stop and sbottom decays at LEP2. The 
dominant decay of the stop and of the sbottom are ti — > eg and bi — > bg [6] when the 
gluino is lighter than the squarks, as in the stable gluino scenario. The branching ratios 
of these decay channels were taken to be 100%. Figure 2 shows the squark production 
and decay diagrams. 

DELPHI data collected in 1994 at a centre-of-mass energy of 91.2 GeV were used to 
search for the process e + e~ — > qqgg. Then, stop and sbottom squarks were searched for 
in DELPHI data collected from 1998 to 2000 at centre-of-mass energies ranging from 189 
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Figure 1: (a) Gluon splitting into a pair of gluinos. (b) Comparison of the cross-section 
(pb) of this process at centre-of-mass energies of 91.2 GeV (LEP1) and of 200 GeV 
(LEP2). 



to 208 GeV. The three analyses presented in this paper are therefore searches for: 



e^ 


-e~ 


-> qqg - 


-> qqgg 


e- 


-e~ 


-> titi - 


c gcg 


e" 


-e" 


— > bibi - 


-> bgbg 



all giving the same topology of two standard jets plus two gluino jets. The gluino could 
either fragment to neutral R° states (gg, guu, ...) or to charged R/ 4 states (gud,...). If 
P is the probability that a gluino fragments to a charged R-hadron, then for P = 1, 
R-hadrons are identified by an anomalous ionizing energy loss in the tracking chambers, 
and for P = 0, the gluino hadronizes into neutral states which reach the calorimeters 
where they deposit a part of their energy. The missing energy carried away by the LSP 
is reduced compared with usual SUSY models with R-parity conservation. For neutral 
R-hadrons, the estimate of the experimental sensitivity depends on the model used to 
calculate the energy loss in the calorimeters. 



2 The DELPHI detector 

The description of the DELPHI detector and its performance can be found in refer- 
ences [7,8]. We only summarize here the parts relevant to the analysis. 

Charged particles are reconstructed in a 1.2 T magnetic field by a system of cylindrical 
tracking detectors. The closest to the beam is the Vertex Detector (VD) which consists 
of three cylindrical layers of silicon detectors at radii 6.3 cm, 9.0 cm and 11.0 cm. They 
measure coordinates in the R<fi plane . In addition, the inner and the outer layers are 
double-sided giving also a z measurement. VD is the barrel part of the Silicon Tracker 

lr The DELPHI coordinate system is defined with z along the e — beam direction; 8 and (f> are the polar and azimuthal 
angles and R is the radial distance from the z axis. 
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Figure 2: Stop (a) and sbottom (b) production and decay at LEP2. 



(ST), which extends the polar angle acceptance down to 10 degrees. The Inner Detector 
(ID) is a drift chamber with inner radius 12 cm and outer radius 22 cm covering polar 
angles between 15° and 165°. The principal tracking detector of DELPHI is the Time 
Projection Chamber (TPC). It is a cylinder of 30 cm inner radius, 122 cm outer radius 
and 2.7 m length. Each end-plate is divided into 6 sectors, with 192 sense wires to al- 
low the dE/dx measurement, and with 16 circular pad rows which provide 3-dimensional 
track reconstruction. The TPC covers polar angles from 20° to 160°. Finally, the Outer 
Detector (OD) consists of drift cells at radii between 192 cm and 208 cm, covering polar 
angles between 43° and 137°. In addition, two planes of drift chambers perpendicular to 
the beam axis (Forward Chambers A and B) are installed in the endcaps covering polar 
angles 11° < 6 < 33° and 147° < 6 < 169°. 

The electromagnetic calorimeters are the High density Projection Chamber (HPC) in 
the barrel region (40° < 9 < 140°) and the Forward Electromagnetic Calorimeter (FEMC) 
in the endcaps (11° < 9 < 36° and 144° < 9 < 169°). In the forward and backward re- 
gions, the Scintillator Tile Calorimeter (STIC) extends the coverage down to 1.66° from 
the beam axis. The number of radiation lengths are respectively 18, 20 and 27 in the 
HPC, the FEMC and the STIC. In the gap between the HPC and the FEMC, hermetic- 
ity taggers made of single layer scintillator-lead counters are used to veto events with 
electromagnetic particles which would otherwise escape detection. Between the HPC 
modules, gaps at 9 — 90° and gaps in are also instrumented with such taggers. Finally, 
the hadron calorimeter (HCAL) covers polar angle between 11° < 9 < 169°. The iron 
thickness in the HCAL is 110 cm which corresponds to 6.6 nuclear interaction lengths. 



3 Data and Monte Carlo samples 

The total integrated luminosity collected by the DELPHI detector in 1994 at the Z° 
peak (y/s = 91.2 GeV) was 46 pb -1 . It corresponded to around 1.6 million hadronic Z° 
events. The Standard Model hadronic background was estimated with the JETSET 7.3 [9] 
program tuned to reproduce LEP1 data [10]. The program described in [1] was used to 
simulate the e + e~ — > qqgg signal. 
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At LEP2, the total integrated luminosity collected by the DELPHI detector at centre- 
of-mass energies from 189 to 208 GeV was 609 pb" 1 . In September 2000, one of the 
twelve TPC sectors (sector 6) stopped functioning. About 60 pb -1 of data were collected 
without that sector until the end of the data taking. The reconstruction programs were 
modified to allow the analysis of the data taken with the DELPHI TPC not fully oper- 
ational. After this modification, there was only a small degradation of the performance 
of the tracking. The data collected in year 2000 were divided in centre-of-mass energy 
windows to optimize the analysis sensitivity. Table 1 summarizes the LEP2 data samples 
used in the analysis. 



Year 


< v/i > (GeV) 


y/s (GeV) 


Integrated luminosity 




Data 


Simulated MC 


(Pb- 1 ) 


1998 


188.6 


189 


158.0 


1999 


191.6 


192 


25.9 




195.6 


196 


76.4 




199.6 


200 


83.4 




201.6 


202 


40.6 


2000 


204.8 


204 


78.1 




206.6 


206 


78.5 




208.1 


208 


7.3 


2000(*) 


206.5 


206.7 


60.6 



Table 1: Total integrated luminosity as a function of the centre-of-mass energy of the 
LEP2 analysed data samples. The third column shows the centre-of-mass energy of the 
simulated events. (*) indicates the data collected by DELPHI in 2000 without the sector 
6 of the TPC. 

The e + e~ interactions leading to four-fermion final states were generated using EX- 
CALIBUR [11]. GRC4F [12] was used to simulate the processes e + e~ — > ez/qq and 
e + e~ — > Z°ee with electrons emitted at polar angles lower than the cut imposed 
in EXCALIBUR. The two-fermion final states were generated with PYTHIA [9] for 
e + e~ — > qq(n7), KORALZ [13] for e + e _ — > t + t~(7), e + e~ — > e + e~ — > 1^(7), 

and BHWIDE [14] for e+e~ -> e+e"^). PYTHIA 6.143 [15] was used to simulate 77 
interactions leading to hadronic final states. BDKRC [16] was used for 77 interactions 
leading to leptonic final states. In all cases, the final hadronization of the particles was 
performed with JETSET [9]. 

The flavour changing decay ti — > eg goes through one-loop diagrams. Therefore, ti 
is expected to be long-lived and to hadronize before its decay. A modified version of 
the SUSYGEN [17] generator was used to simulate this process. Special care was taken 
to introduce hard gluon radiation off the scalar stop at the matrix-element level and 
to treat the stop hadronization as a non-perturbative strong interaction effect. A de- 
tailed description of this hadronization model can be found in [18]. Such a model, based 
on the Peterson function [19], was also used to perform the gluino hadronization into 
R-hadrons. SUSYGEN has also been modified to perform the sbottom decay into bg. 
Figure 3 summarizes the stop and sbottom production and the fragmentation steps. The 
final hadronization was performed using JETSET [9]. 

The Monte Carlo samples used to simulate the Standard Model processes and the 
supersymmetric signals were passed through DELSIM [20], the program simulating the 
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Figure 3: Production and decay of the stop and sbottom squarks. Ellipses indicate the 
color singlet and the color string stretched between the partons. 

full DELPHI detector response. They were subsequently processed with the same recon- 
struction program as the real data. The number of generated events was always several 
times higher than the number expected for the integrated luminosity collected. 



4 R-hadron simulation 



In the analysis, two generic R-hadron states were considered: one charged denoted Pit* 1 
and one neutral, R°, which corresponds to the glueballino, a gg state. It is important to 
understand how an R° would manifest itself in the detector. We refer to the results of 
reference [1]. The energy loss in the scattering on a nucleon R°N — > R°X is given by: 
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where \t\ is the usual momentum transfer invariant for the R° and m x is the mass of the 
system produced in the R°N — > R°X collision. The average energy loss in the reaction 
R°N -> R°X is then given by: 
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From the results of studies in [1], the differential cross-section d ^ mx is taken as: 

da f 1, if \t\ < 1 GeV 2 

d\t\dm x K \ 0, if \t\ > 1 GeV 2 

The average number of collisions of an R° particle in the calorimeters is given by the 
depth of the calorimeter in units of equivalent iron interaction lengths, Ay. In DELPHI, 
the electromagnetic calorimeter's thickness represents around 1 At while this value is 
6.6 Ay for the hadronic calorimeter. We have adopted a correction factor for the inter- 
action length of 9/16 as suggested in reference [1]: a factor Ca/Cf = 9/4 comes from 
the colour octet nature of the R° constituents increasing the aRo N cross-section as com- 
pared to a^N, while a factor {r 2 RO ) / '(r 2 ) = 1/4 takes into account the relative size of the 
R-hadrons as compared to standard hadrons. On average, neutral R-hadrons should un- 
dergo 4.3 collisions in DELPHI calorimeters. Figure 4 shows the total energy loss by an 
R° after 4 collisions in iron. 

The difficulty in separating the signal of a neutral R-hadron from the background 
increases with the amount of energy lost in the calorimeters. The choice of the interac- 
tion model made here is conservative in this respect. The R° scatters were subsequently 
treated in the DELPHI detector simulation as with the energy that the R° should 
deposit in four collisions according to the above formula. 

The charged R-hadrons were treated as heavy muons to reproduce the anomalous 
dE/dx signature. In this case, only the tracking information was used to calculate the 
R-hadron momentum. 




Figure 4: Average energy loss by neutral R-hadrons in the DELPHI calorimeters as a 
function of their initial energy for different mass cases. 
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5 Particle Identification and analysis method 

5.1 Particle Identification and event preselection 

Particle selection was identical for LEP1 and LEP2 data. Reconstructed charged 
particles were required to have momenta above 100 MeV/c with Ap/p < 1, where Ap is 
the momentum error, and impact parameter below 5 cm in the transverse plane and below 
10 cm/ sin 6* in the beam direction. More stringent cuts were applied for tracks without 
TPC information. A cluster in the calorimeters was selected as a neutral particle if not 
associated to a charged particle and if the cluster energy was greater than 500 MeV in the 
HPC, 400 MeV in the FEMC, 300 MeV in the STIC or 900 MeV in the HAC. Particles 
were then clustered into jets with the DURHAM algorithm [21]. b-quarks were tagged 
using a probabilistic method based on the impact parameters of tracks with respect to 
the main vertex. A combined b-tagging variable was defined by including the properties 
of secondary vertices [22]. 

Events were then kept if there were at least two charged particles, and at least one 
with a transverse momentum above 1.5 GeV/c, and if the transverse energy 2 exceeded 
4 GeV. 

5.2 Neural networks 

A neural network allows one discriminating variable to be constructed from the set of 
variables given as input. The form used here contains three layers of nodes: the input 
layer where each neuron corresponds to a discriminating variable, the hidden layer, and 
the output layer which is the response of the neural network. The layers were connected 
in a "feed forward" architecture. The back-propagation algorithm was used to train 
the network with simulated events. This entails minimising a x 2 to adjust the neurons' 
weights and connections. An independent validation sample was also used not to overtrain 
the network. The outputs of neural networks were used to isolate events containing two 
neutral R-hadrons at LEP2. 

5.3 General description of the analyses 

The LEP1 and LEP2 data analyses have all the same final states which consists of 
two jets and two gluino jets. Depending on the probability P that the gluino hadronizes 
into a charged R-hadron, three topologies are possible: 

• charged: e + e" — > qqR ± R ± 

• mixed: e + e~ — > qqR ± R° 

• neutral: e + e~ — > qqR°R° 

For LEP1 data, the search analyses corresponding to the charged and the mixed topolo- 
gies were identical, while they were optimised separately for LEP2 data. For these two 
topologies, the LEP1 and LEP2 analyses were only based on anomalous ionizing energy 
loss, and no quark tagging was applied. 

For the neutral topology, LEP1 data were analyzed using sequential cuts without try- 
ing to tag the quark flavor. Since there are b-quarks in the final state of the sbottom 
decay, the b-quark tagging improves the isolation of this signal from the Standard Model 
backgrounds. Therefore, the b-tagging variable was used in the stop and sbottom search 

2 The transverse energy is denned as the sum of J p\ + m? over all particles; pj< is the transverse momentum. 
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analyses at LEP2 in the neutral topology. This variable was included into the neural 
networks which were used to isolate the stop and the sbottom signals. 

6 Search for a stable gluino at LEP1 

6.1 Search for qqR^Rr 11 and qqP^R events 

The same analysis based on the dE/dx measurement was performed to identify 
qqR ± R ± and qqT^R events. 

In the preselection step, events were required to contain at least 5 charged particles. 
At least one of these had to satisfy the following conditions. The track was required to be 
reconstructed including a TPC track element and to have a momentum above 10 GeV/c. 
At least 80 wires of the TPC were required to have been included in the dE/dx mea- 
surement. The dE/dx had to be either greater than 1.8 mip (units of energy loss for 
a minimum ionizing particle), or less than the dE/dx expected for a particle of mass 
equal to 1 GeV/c 2 . The V 23 variable is the y cut value in the DURHAM algorithm for 
which the number of jets changes between two and three. qqR ± R ± and qqR ± R° events 
contain three or four jets. Thus, F 23 was required to be less than 0.01. Figure 5 shows a 
comparison between simulated and real data at this level. 

A R^ candidate had to satisfy the following conditions: it had to be reconstructed 
with the VD, the ID and the TPC detectors, and the dE/dx measurement had to be 
based on at least 80 wires of the TPC. In addition, the energy of the other particles in 
a 15° cone around the R^ candidate had to be less than 2 GeV. Finally, its associated 
electromagnetic energy had to be less than 5 GeV. 

The final selection was performed by cuts in the plane (P,dE/dx). Figure 6 shows the 
expected dE/dx as a function of the particle momentum. The analysis was separated 
into two mass windows: 

• mg < 14 GeV/c 2 : 

Here, charged R-hadrons were identified by low dE/dx values. The R^ candidates 
were selected if their momentum was greater than 15 GeV/c, and if their dE/dx was 
less than the dE/dx expected for a particle of mass equal to 3 GeV/c 2 . 

• mg > 14 GeV/c 2 : 

In this mass window, R-hadrons were identified by high dE/dx values. The R^ 1 
candidates were selected if their dE/dx was greater than 2 mip. 

The final selection was performed by requiring at least one charged R-hadron candidate 
in either mass window. Table 2 contains the number of events selected after each cut 
of this analysis. For mg < 14 GeV/c 2 , 5 events were selected when 4.2 were expected. 
These numbers are 12 and 13.5 in the mg > 14 GeV/c 2 mass window. Unlike the expected 
signal, all selected candidates in the data have only one particle with anomalous dE/dx. 
Figure 7 shows the signal detection efficiencies. For qqR ± R ± , they ranged from a few 
percent for gluino masses close to 2 GeV/c 2 to around 50% for gluino masses of the order 
of 25 GeV/c 2 . qqR ± R° efficiencies were about half of the qqR ± R ± ones. 

6.2 Search for qqR°R° events 

The search for qqR°R° events was performed at LEP1 with a sequential cut analysis. 
It was based on the search for the small part of missing energy carried away by the neutral 
R-hadrons. Hadronic events were first selected by requiring at least 5 charged particles. 
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After forcing the events into two jets, the acollinearity 3 was required to be greater than 
20° to reduce the huge number of background Standard Model Z° — > qq events. 

The following cuts were applied to reduce the number of hadronic 77 interactions. 
The number of tracks reconstructed with the TPC had to be greater than 4, and the 
energy of the particles with tracks reconstructed using only the VD and ID detectors had 
to be less than 20% of the total energy. The energies in 40° and 20° cones around the 
beam axis were required to be less than 40% and 10% of the total energy respectively. 
The transverse energy had to be greater than 20 GeV. 

Hadronic events with missing energy were then selected in the barrel region of the 
detector. The visible mass was required to be less than 60 GeV/c 2 . The thrust axis 
and the missing momentum had to point in the polar regions [37°, 143°] and [45°, 135°] 
respectively. Figure 8 shows a comparison between data and simulation at this level of 
the selection. 

The Y23 quantity was then required to be less than 0.01 and events had to contain less 
than 20 charged particles. In order to reduce the number of events with two back-to-back 
jets, the acoplanarity 4 was required to be greater than 10° and the thrust to be less than 
0.95. The final cut was bi-dimensional. The value of the variable Mj e fi/ Ej e n-\- Mj^/ -Ejet2 
was calculated from the two jets reconstructed with the DURHAM algorithm. Events 
were rejected if this variable was greater than 0.45 and if the acollinearity was less than 
50°. Table 3 shows the number of events after each cut of the qqR°R° analysis. 12 events 
were selected in the data while 10.6 were expected in the hadronic background. Signal 
efficiencies as a function of the gluino mass are shown in figure 7. They ranged from a 
few percent for low gluino masses to around 20% for rrig = 18 GeV/c 2 . 

7 Search for a stable gluino at LEP2 

7.1 Preselection 

A common preselection for the charged and neutral R-hadron analyses was applied to 
reduce the background coming from soft 77 interactions. The cuts are the same for the 
stop and sbottom analysis at all centre-of-mass energies ranging from 189 to 208 GeV. 

First, events were forced into two jets. To select hadronic events, the number of 
charged particles reconstructed with the TPC was required to be greater than three, and 
the energy in the STIC to be less than 70% of the total visible energy. The polar angle of 
the thrust axis had to be in the interval [20°, 160°]. Then, quality cuts were applied. The 
fraction of good tracks was defined as the ratio between the number of charged particles 
remaining after the track selection divided by this number before the selection. It had 
to be greater than 35%. In addition, the scalar sum of particle momenta reconstructed 
with the TPC was required to be greater than 55% of the total reconstructed energy 
and the number of charged particles to be greater than six. To remove radiative events, 
the energy of the most energetic neutral particle had to be less than 40 GeV. Table 4 
contains the number of events after each of these cuts. 

For the qqR ± R ± and qqR ± R° analyses, charged R-hadron candidates were defined at 
this level. They had to be reconstructed with the VD, ID and TPC detectors and their 
momentum was required to be greater than 10 GeV/c. At least 80 sense wires of the TPC 
were required to have contributed to the measurement of their dE/dx. Their associated 
electromagnetic energy was required to be less than 5 GeV, and the energy of the other 

3 The acollinearity of two jets is defined as the complement of the angle between their directions. 

4 The acoplanarity of two jets is defined as the complement of the angle between their directions projected onto the 
plane perpendicular to z. 
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charged particles in a 15° cone around a candidate had to be less than 5 GeV. In 2000 , 
the dE/dx could not be used in sector 6 of the TPC for almost any of the data. For this 
sample, charged R-hadron candidates in this sector were removed. 

7.2 Search for qqR^Rr 11 events 

The search for qqR ± R ± events was exactly the same for the stop and the sbottom 
analyses. Events were selected if they contained at least two charged R-hadron candi- 
dates. Figure 9 shows the momentum and the dE/dx distribution of the selected 
candidates. Table 4 shows the number of selected events. 

The analysis was then separated into three windows in gluino mass, and cuts in the 
plane (P, dE/dx) were applied: 

• m g < 30 GeV/c 2 : 

events had to contain at least one charged R-hadron candidate with momentum 
greater than 20 GeV/c, and with dE/dx less than the dE/dx expected for a particle 
of mass equal to 3 GeV/c 2 . 

• 30 GeV/c 2 < m g < 60 GeV/c 2 : 

events were selected if they contained at least two charged R-hadron candidates 
with dE/dx both greater than the dE/dx expected for a particle of mass equal 
to 30 GeV/c 2 , and less than the dE/dx expected for a particle of mass equal to 
60 GeV/c 2 . Moreover, this dE/dx had also to be either less than the dE/dx expected 
for a particle of mass equal to 1 GeV/c 2 , or greater than 1.8 mip. 

• m g > 60 GeV/c 2 : 

events were kept if they contained at least two charged R-hadron candidates with 
dE/dx greater than the dE/dx expected for a particle of mass equal to 60 GeV/c 2 . 

In all LEP2 data which were analysed, no events were selected in any of these windows. 
The number of expected Standard Model background events were 0.115, 0.009 and 0.011 
in the analyses for m g < 30 GeV/c 2 , 30 GeV/c 2 < m g < 60 GeV/c 2 and m g > 60 GeV/c 2 
respectively. Table 5 contains the number of events expected for the different centre-of- 
mass energies. Figure 10 shows the signal detection efficiencies near the kinematical limit 
(m^ = 90 GeV/c 2 ). The difference between stop and sbottom efficiencies is not large. 
The highest efficiencies were always obtained for high gluino masses, where the dE/dx is 
very high. 

7.3 Search for qqR^R events 

The search for qqR ± R° events was also the same for the stop and sbottom analyses. 
Events were selected if they contained at least one charged R-hadron candidate. Figure 11 
shows the momentum and the dE/dx distribution of the selected candidates. Table 4 
shows the number of selected events. 

The analysis was then separated into three gluino mass windows, and cuts in the plane 
(P,dE/dx) were applied: 

• m t < 30 GeV/c 2 : 

events had to contain at least one charged R-hadron candidate with momentum 
greater than 20 GeV/c, and with dE/dx less than the dE/dx expected for a particle 
of mass equal to 3 GeV/c 2 . 

• mg > 60 GeV/c 2 : 

events were kept if they contained at least one charged R-hadron candidate with 
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dE/dx greater than the dE/dx expected for a particle of mass equal to 60 GeV/c 2 , 
and greater than 2 mip. 

• 30 GeV/c 2 < wig < 60 GeV/c 2 : 

events selected in either of the above windows (higher or lower trig) were accepted. 

Three, nine and six events were selected in the mass windows m g < 30 GeV/c 2 , 
30 GeV/c 2 < rag < 60 GeV/c 2 and m g > 60 GeV/c 2 respectively. The number of 
expected background events were 1.6, 8.2 and 6.6. All selected events in the data are 
more likely Standard Model instead than signal like. In particular, they do not follow any 
mass iso-curve in the (P, dE/dx) plane. Table 6 contains the number of selected events as 
a function of the centre-of-mass energy. Figure 12 shows the signal detection efficiencies 
near the kinematical limit (m^ = 90 GeV/c 2 ). The highest efficiencies were obtained for 
high gluino masses where the dE/dx is very high. 

7.4 Search for qqR°R° events 

After the preselection described in section 7.1, the transverse missing momentum was 
required to be greater than 4 GeV/c, the angle of the missing momentum had to point in 
the polar angle region [20°, 160°], and the energy in a 40° cone around the beam axis was 
required to be less than 40% of the event energy. A veto algorithm was then applied based 
on the hermeticity taggers at polar angles close to 40° and 90°. Figures 13 and 14 show 
data Monte Carlo comparisons following this selection and table 7 gives the observed and 
expected event numbers at the different steps. 

The stop and sbottom analyses were then separated for different ranges of the mass 
difference Am between the squark and the gluino: 

• Am < 20 GeV/c 2 : 

For high gluino masses, the energy deposited by the neutral R-hadrons is quite small. 
In this respect, the gluino is not so different from a neutralino, and the qi — > qg 
events resemble qi — > qx? events. 

• Am > 20 GeV/c 2 : 

In this case, the gluino deposits more energy. 

The neural networks were trained to isolate the qqR°R° signal in both Am windows. 
They were trained separately on stop signals or sbottom signals. The neural network 
structure was the same for the stop and the sbottom searches. It consisted of 10 input 
nodes, 10 hidden nodes and 3 output nodes. 

For Am < 20 GeV/c 2 , the 10 input variables were: the ratio between the transverse 
missing momentum and the visible energy, the transverse energy, the visible mass, the 
softness defined as Mj et \/ Ej et i+Mj et2 / Ej et2 , the acollinearity, the quadratic sum of trans- 
verse momenta of the jets \J (P/ e<1 ) 2 + (P/ e<2 ) 2 , the acoplanarity, the sum of the first and 
third Fox- Wolfram moments, the polar angle of the missing momentum and finally the 
combined b-tagging probability [22]. 

For Am > 20 GeV/c 2 , the 10 input variables were: the charged energy, the transverse 
charged energy, the visible mass, the thrust, the effective centre-of-mass energy [23], the 
acollinearity, the acoplanarity, the sum of the first and third Fox- Wolfram moments, the 
sum of the second and fourth Fox- Wolfram moments, and finally the combined b-tagging 
probability. 

The neural networks were trained to discriminate the signal from the combined two- 
fermion and four-fermion backgrounds, and from the 77 interactions leading to hadronic 
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final states. The first output node was trained to identify both Z — > qq and four-fermion 
events. The second node identified the 77 interactions leading to hadronic final states. 
And the third node was trained to select the signal. The three output nodes were useful 
in the training of the network, but the selection was made according to the output of 
the signal node only. Figure 15 shows the number of events as a function of the signal 
efficiency for the two mass analysis windows of the stop and the sbottom analysis. The 
number of real events was in agreement with the Standard Model predictions over the full 
range of the neural network outputs. The optimisations of the final cuts were performed 
by minimising the expected confidence level of the signal hypothesis [24]. 

Tables 8 and 9 contain the numbers of events selected in the stop and sbottom anal- 
yses. Combining all data from 189 to 208 GeV, 32 and 11 events were selected in the 
stop analysis for Am > 20 GeV/c 2 and Am < 20 GeV/c 2 , while the expected number 
of events were 30.1 and 11.1. In the sbottom analysis, no candidates were observed for 
Am > 20 GeV/c 2 and five were selected for Am < 20 GeV/c 2 . The expected number of 
events were 3.0 and 5.3. Figure 16 shows the signal detection efficiencies for the stop and 
for the sbottom (m^ = 90 GeV/c 2 ). They are very low when the gluino mass is close 
to zero. 



8 Results 

No excess of events was observed in any analysis performed at LEP1 or at LEP2 in the 
stable gluino scenario. Results were therefore combined to obtain excluded regions at 95% 
confidence level in the parameter space. The limits were computed using the likelihood 
ratio method described in [24]. For different values of the parameter P describing the 
probability that the gluino hadronizes to a charged R-hadron, the relative cross-sections 
for the different channels were given by: 

f (j(R ± R ± ) = P 2 a 

I ( j(R ± R°) = 2P(1-P)a 

( a(R°R°) = (l-P)V 

where a was either the e + e" — > qqgg cross section at LEP1, the titi or the bibi cross- 
section at LEP2. 

For the LEP1 analysis, results were interpreted in terms of excluded gluino masses for 
different P. Figure 17 shows the excluded region at 95% confidence level. From this figure, 
a stable gluino with mass between 2 and 18 GeV/c 2 is excluded regardless of the charge 
of the R-hadrons. The minimum upper limit, 18 GeV/c 2 , is obtained for intermediate 
values of P (between 0.2 and 0.45), while an upper limit of 23 GeV/c 2 is obtained for 
P=l. 

For the LEP2 analysis, excluded regions in the planes {m- tl ,mg) and (m^,mg) were 
derived for five different values of P: 0., 0.25, 0.5, 0.75 and 1. Moreover, the stop and 
sbottom cross-sections were calculated for two cases. In the first case, the squark mixing 
angle was set to zero which corresponds to the maximal cross-sections. In the second 
case, the mixing angle was equal to 56° for the stop and to 68° for the sbottom, which 
gives the minimal cross-sections. Squark masses below 50 GeV/c 2 were not taken into 
account in these analyses. Figures 18 and 19 show the excluded regions thus obtained. 
In the case of minimal cross-sections, a hole appears in the sbottom exclusion histograms 
around (60 GeV/c 2 , 50GeV/c 2 ) in the plane (m^^mg) for intermediate values of P. It can 
be explained by different effects: 



13 



• In the analyses searching for charged R-hadrons, the signal detection efficiencies are 
small in this region of the (mq^wig) plane, because the dE/dx distribution of the 
charged R-hadrons as a function of the momentum crosses the band corresponding 
to Standard Model particles (cf Figure 6). 

• When the Z°— qj — q 1 coupling is suppressed, the squark pair production cross-section 
is much more reduced for the sbottom than for the stop. 

• The visible energy of qqR°R° events becomes small for low values of AM. At 
AM = 5 GeV/c 2 , this has no effect on the stop exclusion results, because of the 
c-quark which is lighter than the b-quark. 

Lower limits on the stop and sbottom masses are given in table 10 for Am > 10 GeV/c 2 
and for a gluino mass greater than 2 GeV/c 2 for different values of P. 

9 Conclusion 

The analysis of the LEP1 data collected in 1994 excludes at 95% confidence level a 
stable gluino with mass between 2 and 18 GeV/c 2 . These limits are valid for any charge 
of the produced R-hadrons. 

Stop and sbottom squarks have been searched for in the 609 pb -1 collected by DELPHI 
at centre-of-mass energies ranging from 189 to 208 GeV. In the stable gluino scenario, 
the dominant decays are ti — > eg and b x — > bg. No deviation from Standard Model 
predictions was observed. The observed limits at 95% confidence level are: 

• mi Y > 90 GeV/c 2 , and rag > 96 GeV/c 2 for purely left squarks. 

• m- tl > 87 GeV/c 2 , and > 82 GeV/c 2 independent of the mixing angle. 
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Data 


Z u — > qq 
background 


anomalous dE/dx 
1 candidate 


99322 
24566 
421 


97170 ± 200 
25794 ± 104 
464 ± 14 


low dE/dx (wig< 14) 
high dE/dx (m t > 14) 


5 

12 


4.2 ± 1.3 
13.5 ± 2.4 



Table 2: Number of events selected after each cut of the charged R-hadron analysis at 
LEP1. 





Data 


Z u — > qq 
background 


Acolinearity 


41231 


34853 ± 120 


Ntpc 


38977 


33807 ± 120 


%Evd~id 


36877 


32419 ± 120 


E40/ E vis 


19309 


15311 ± 80 




16664 


13480 ± 75 


E t 


16317 


13453 ± 75 


M vis 


5932 


6353 ± 52 


COS 9 thrust 


5384 


5725 ± 49 


cos 6 Pmis \ 


2527 


2294 ± 31 




214 


194 ± 9 


■Nchar. 


183 


161 ± 8 


Acoplanarity 


134 


115 ± 7 


Thrust 


105 


81.7 ± 5.9 


Acol. vs M jetl /E jetl + M jet2 /E jet2 


12 


10.6 ± 2.1 



Table 3: Number of events selected after each cut of the qqR°R° analysis at LEP1. 



1G 



Cuts 



Data 



Simulation 



4-fermions 



2-fermions 



77 



Ntpc 
Estic/ E v 

^thrust 

Ptpc 
N cha . 

rpmax. 
neu. 



175436 
145810 
54838 
48475 
45816 
41880 



164146 ± 
141362 ± 
54933 ± 
48846 ± 
46227 ± 
42113 ± 



105 
95 
45 
43 
37 
37 



1 candidate 



12418 


± 


13 


12062 


± 


13 


9739 


± 


10 


9141 


± 


10 


9040 


± 


9 


8802 


± 


9 



50391 ± 
48170 ± 
31510 ± 
27580 ± 
26969 ± 
23108 ± 



24 
24 
23 
23 
16 
16 



101338 ± 102 



81131 ± 
13685 ± 
12126 ± 
10219 ± 
10203 ± 



91 
38 
35 
32 
32 



2187 



2279.1 ± 6. 



1746.8 ± 4.6 



470.3 ± 3.4 



62.6 ± 2.6 



2 candidates 



74 



79.2 ± 0.8 



75.2 ± 0.7 



3.4 ± 0.3 



0.5 ± 0.2 



Table 4: Number of events after each cut of the LEP2 preselection. 189 to 208 GeV data 
are added. 





wig < 30 GeV/c 2 


30 GeV/c 2 < m g < 60 GeV/c 2 


wig > 60 GeV/c 2 




Data 


Simulation 


Data 


Simulation 


Data 


Simulation 


188.7 





0.029±0.016 





0.001±0.001 





0.003±0.003 


191.6 





0.005±0.004 





0.001±0.001 





0.001±0.001 


195.6 





0.025±0.020 





0.001±0.001 





0.001±0.001 


199.6 





0.007±0.007 





0.001±0.001 





0.001±0.001 


201.7 





0.011±0.007 





0.001±0.001 





0.001±0.001 


204.8 





0.009±0.009 





0.001±0.001 





0.001±0.001 


206.7 





0.012±0.009 





0.001±0.001 





0.001±0.001 


208.1 





0.000±0.000 





0.001±0.001 





0.001±0.001 


206.5(*) 





0.017±0.013 





0.001±0.001 





0.001±0.001 


total 





0.115±0.033 





0.009±0.003 





0.011±0.004 



Table 5: Number of events selected by the qq^R 1 * 1 analysis at LEP2. (*) indicates 2000 
data taken without sector 6 working. 





m g < 30 GeV/c 2 


30 GeV/c 2 < m g < 60 GeV/c 2 


m g > 60 GeV/c 2 




Data 


Simulation 


Data 


Simulation 


Data 


Simulation 


188.7 





0.577±0.101 





2.514±0.264 





1.937±0.243 


191.6 





0.030±0.011 


1 


0.454±0.132 


1 


0.425±0.131 


195.6 


2 


0.135±0.042 


2 


0.783±0.097 





0.648±0.088 


199.6 





0.266±0.071 


1 


1.333±0.158 


1 


1.068±0.141 


201.7 





0.097±0.025 


2 


0.489±0.056 


2 


0.392±0.050 


204.8 


1 


0.208±0.051 


1 


0.961±0.106 





0.753±0.093 


206.6 





0.187±0.040 


1 


0.848±0.089 


1 


0.661±0.079 


208.1 





0.011±0.005 





0.095±0.014 





0.085±0.013 


206.5(*) 





0.083±0.025 


1 


0.679±0.074 


1 


0.596±0.070 


total 


3 


1.59±0.15 


9 


8.16±0.39 


6 


6.57±0.36 



Table 6: Number of events selected by the qql^R analysis at LEP2. (*) indicates 2000 
data taken without sector 6 working. 
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Cuts 


Data 


simulation 


4-fermions 


2-fermions 


77 


pmis. 

o P ■ 

1 mis. 

-E40 / E V i S , 
Hermeticity 


26423 
16379 
14694 
14422 


26938 ± 20 
16821 ± 15 
15231 ± 14 
14651 ± 14 


8117 ± 8 
7191 ± 6 
6395 ± 6 
6150 ± 6 


18012 ± 15 
9088 ± 12 
8471 ± 12 
8140 ± 12 


809 ± 9 
542 ± 8 
364 ± 6 
361 ± 6 



Table 7: qqR°R° analysis at LEP2: number of events after each selection cut. Data with 
y/s in the range 189 GeV-208 GeV are included. 





AM > 20 GeV/c 2 


AM < 20 GeV/c 2 




Data 


Simulation 


Data 


Simulation 


188.7 


4 


6.634±0.741 


6 


3.685±1.158 


191.6 


4 


1.054±0.115 





0.482±0.097 


195.6 


5 


3.532±0.236 


3 


1.408±0.256 


199.6 


7 


4.324±0.270 





1.617±0.280 


201.7 


1 


2.055±0.140 





0.836±0.138 


204.8 


4 


4.432±0.302 





1.197±0.268 


206.6 


5 


4.287±0.290 


2 


1.227±0.272 


208.1 





0.418±0.031 





0.117±0.026 


206.5(*) 


2 


3.411±0.203 





0.556±0.072 


total 


32 


30.2±1.0 


11 


11.1±1.3 



Table 8: Number of events selected by the stop qqR°R° analysis. (*) indicates 2000 data 
taken without sector 6 working. 





AM > 20 GeV/c 2 


AM < 20 GeV/c 2 




Data 


Simulation 


Data 


Simulation 


188.7 





1.038±0.657 


3 


2.475±1.139 


191.6 





0.085±0.026 





0.180±0.058 


195.6 





0.314±0.088 


1 


0.571±0.178 


199.6 





0.295±0.094 





0.663±0.195 


201.7 





0.225±0.049 


1 


0.319±0.096 


204.8 





0.295±0.040 





0.428±0.174 


206.6 





0.312±0.05 





0.379±0.175 


208.1 





0.022±0.005 





0.037±0.017 


206.5(*) 





0.417±0.054 





0.281±0.160 


total 





3.0±0.7 


5 


5.3±1.2 



Table 9: Number of events selected by the sbottom qqR°R° analysis. (*) indicates 2000 
data taken without sector 6 working. 
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Stop 


Sbottom 


p 


0i = o° 


Of 


= 56° 


% = o° 


9~ b = 68° 


0.00 


92 




87 


98 


86 


0.25 


90 




87 


96 


82 


0.50 


92 




89 


94 


82 


0.75 


94 




92 


94 


84 


1.00 


95 




94 


95 


87 



Table 10: Upper limits on the stop and sbottom masses as a function of the probability 
P that the gluino hadronizes to charged R-hadrons. Limits are set for Am > 10 GeV/c 2 
and for a gluino mass greater than 2 GeV/c 2 . Mixing angles equal to zero corresponds 
to purely left-handed squarks, while 6> t ~ = 56° and 9^ = 68° corresponds to minimal 
cross-section case. 
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Figure 5: Comparison between data and simulation in the search for charged R-hadrons 
at LEP1. The plots show characteristic distributions before the selection of the charged 
R-hadron candidates: (a) the momentum, (b) the total energy in a 15° degree cone 
around the particle, (c) its electromagnetic and (d) hadronic energy. Dotted lines show the 
qqR ± R ± and qqR ± R° signal distributions with arbitrary normalization when all simulated 
samples are added together. 
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10 20 30 40 50 60 70 80 90 100 

Momentum (GeV/c) 



Figure 6: The lines show the theoritical ionization energy loss as a function of the momen- 
tum of the particle for different mass cases (GeV/c 2 ). The black points are reconstructed 
tracks coming from Z° — > qq events generated at y/s =200 GeV, while the stars and 
the triangles correspond to R 1 * 1 tracks with mass of 10 and 70 GeV/c 2 respectively. The 
dE/dx is expressed in units of energy loss for a minimum ionizing particle. 



21 



O 

g 0.7 
o 



0.6 



0.5 



0.4 



0.3 



DELPHI: LEP1 






■■■■ qqR +/ "R" /+ signal 




* 

• 

.« 


qqR R signal 
— qqR°R° signal 




* 

• 

* 

* 

• 

* 

* 

• 




• 
* 


* 






* 

* 

■ 

• 

* 

* 

* 

■ 




* ♦ 


* 

• 
* 






Gluino mass (GeV/c ) 



Figure 7: Signal detection efficiencies in the LEP1 data analyses as a function of the 
gluino mass for the three signal topologies qqR R , qqR ± R° and qqR°R°. The search 
for charged R-hadrons was optimized separately for low gluino masses (< 14 GeV/c 2 ) 
and for high gluino masses (> 14 GeV/c 2 ). 
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Figure 8: Comparison between data and simulation in the qqR°R° analysis at LEP1. (a) 
visible energy, (b) visible mass, (c) acoplanarity (d) DURHAM distance Y23. Dotted lines 
show the qqR°R° signal distributions with arbitrary normalization when all simulated 
samples are added together. 
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DELPHI: 189 to 208 GeV data 
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Figure 9: Momentum and dE/dx of the charged R-hadron candidates selected by the 
qqR ± R ± analysis at LEP2. Data taken in the centre-of-mass energy range between 189 
and 208 GeV were included. Right-hand side histograms show the expected distributions 



for one possible stop and sbottom signal (m. 
y/s = 200 GeV with arbitrary normalization. 
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Figure 10: Signal detection efficiencies at =200 GeV for the stop (a) and sbottom 
(b) qqR ± R ± analysis as a function of the gluino mass (m^ = 90 GeV/c 2 ). AM is the 
mass difference between the squark and the gluino. Vertical lines show the limits of the 
mass analysis window. The last one ends at AM = 5 GeV/c 2 which corresponds the last 
simulated signal points. 
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Figure 11: Momentum and dE/dx of the charged R-hadron candidates selected by the 
qql^R analysis at LEP2. Data taken in the centre-of-mass energy range between 189 
and 208 GeV were included. Right-hand side histograms show the expected distributions 



for one possible stop and sbottom signal (m. 
yfs = 200 GeV with arbitrary normalization. 
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Figure 12: Signal detection efficiencies at -Js =200 GeV for the stop (a) and sbottom 
(b) qqR ± R° analysis as a function of the gluino mass (m^ = 90 GeV/c 2 ).AM is the 
mass difference between the squark and the gluino. Vertical lines show the limits of the 
mass analysis window. The last one ends at AM = 5 GeV/c 2 which corresponds the last 
simulated signal points. 
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Number of Particles Number of Particles 

Figure 13: Data-simulation comparison at the preselection level of the LEP2 qqR°R° 
analysis. Data taken in the centre-of-mass energy range between 189 and 208 GeV were 
included. Right-hand side histograms show the expected distributions with arbitrary 
normalization for the stop and the sbottom signal at 200 GeV when all simulated samples 
are added together. 
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Figure 14: Data-simulation comparison at the preselection level of the LEP2 qqR°R° 
analysis. Data taken in the centre-of-mass energy range between 189 and 208 GeV were 
included. Right-hand side histograms show the expected distributions with arbitrary 
normalization for the stop and the sbottom signal at 200 GeV when all simulated samples 
are added together. 
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Figure 15: Numbers of events as a function of the signal efficiencies for the stop and sbot- 
tom analysis. Data taken in the centre-of-mass energy range between 189 and 208 GeV 
were included, (a) stop analysis for Am > 20 GeV/c 2 and (b) for Am < 20 GeV/c 2 , (c) 
sbottom analysis for Am > 20 GeV/c 2 and (b) for Am < 20 GeV/c 2 . The squark and 
gluino mass values used for the signal detection efficiencies are indicated on the x axis, 
(mq,mg). 
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Figure 16: Signal detection efficiencies at yfs =200 GeV for the stop (a) and sbottom (b) 
qqR°R° analysis as a function of the gluino mass (m^ = 90 GeV/c 2 ). 
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Figure 17: Results of the LEP1 analysis: excluded region at 95% confidence level in the 
plane (mg,P). P is the probability that the gluino hadronizes to a charged R-hadron. The 
shaded region corresponds to the observed exclusion and the line to the expected one. 
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Figure 18: Results of the LEP2+LEP1 stop analysis: excluded region at 95% confidence 
level in the plane (m^mg). The line corresponds to the exclusion for purely left stop, 
and the shaded region to exclusion obtained for the mixing angle giving the minimal 
cross-section. Excluded regions are given for different values of P, the probability that 
the gluino hadronizes to charged R-hadron: 0, 0.25, 0.5, 0.75 and 1. 
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Figure 19: Results of the LEP2+LEP1 sbottom analysis: excluded region at 95% con- 
fidence level in the plane (mg ,m.g). The line corresponds to the exclusion for purely 
left sbottom, and the shaded region to exclusion obtained for the mixing angle giving the 
minimal cross-section. Excluded regions are given for different values of P, the probability 
that the gluino hadronizes to charged R-hadron: 0, 0.25, 0.5, 0.75 and 1. 



